source('~/Dropbox/cloud_work/papers in progress/death_free_market/Supplementary Material/Figures/manager_energy_plot.R')
source('~/Dropbox/cloud_work/papers in progress/death_free_market/Supplementary Material/Figures/manager_energy_plot.R')
source('~/Dropbox/cloud_work/papers in progress/death_free_market/Supplementary Material/Figures/manager_energy_plot.R')
source('~/Dropbox/cloud_work/papers in progress/death_free_market/Supplementary Material/Figures/manager_energy_plot.R')
source('~/Dropbox/cloud_work/papers in progress/death_free_market/Supplementary Material/Figures/landscape_inset.R')
source('~/Dropbox/cloud_work/papers in progress/death_free_market/Supplementary Material/Figures/manager_energy_model_best.R')
source('~/Dropbox/cloud_work/papers in progress/death_free_market/Supplementary Material/Figures/government_energy_model_rainbow.R')
source('~/Dropbox/cloud_work/papers in progress/death_free_market/Supplementary Material/Figures/government_energy_model_rainbow.R')
source('~/Dropbox/cloud_work/papers in progress/death_free_market/Supplementary Material/Figures/government_energy_model_rainbow.R')
source('~/Dropbox/cloud_work/papers in progress/death_free_market/Supplementary Material/Figures/government_energy_model_rainbow.R')
source('~/Dropbox/cloud_work/papers in progress/death_free_market/Supplementary Material/Figures/government_energy_model_rainbow.R')
source('~/Dropbox/cloud_work/papers in progress/death_free_market/Supplementary Material/Figures/government_energy_model_rainbow.R')
source('~/Dropbox/cloud_work/papers in progress/death_free_market/Supplementary Material/Figures/government_energy_model_rainbow.R')
source('~/Dropbox/cloud_work/papers in progress/death_free_market/Supplementary Material/Figures/government_energy_model_best.R')
source('~/Dropbox/cloud_work/papers in progress/death_free_market/Supplementary Material/Figures/RUN_ALL_FIGURES.R')
source('~/Dropbox/cloud_work/papers in progress/death_free_market/Supplementary Material/Figures/manager_energy_model_rainbow.R')
source('~/Dropbox/cloud_work/papers in progress/death_free_market/Supplementary Material/Figures/energy_government_model.R')
source('~/Dropbox/cloud_work/papers in progress/death_free_market/Supplementary Material/Figures/energy_government_model.R')
View(sim_best)
source('~/Dropbox/cloud_work/papers in progress/death_free_market/Supplementary Material/Figures/energy_government_model.R')
source('~/Dropbox/cloud_work/papers in progress/death_free_market/Supplementary Material/Figures/energy_government_model.R')
source('~/Dropbox/cloud_work/papers in progress/death_free_market/Supplementary Material/Figures/energy_government_model.R')
source('~/Dropbox/cloud_work/papers in progress/death_free_market/Supplementary Material/Figures/energy_government_model.R')
source('~/Dropbox/cloud_work/papers in progress/death_free_market/Supplementary Material/Figures/energy_government_model.R')
source('~/Dropbox/cloud_work/papers in progress/death_free_market/Supplementary Material/Figures/energy_government_model.R')
source('~/Dropbox/cloud_work/papers in progress/death_free_market/Supplementary Material/Figures/energy_government_model.R')
source('~/Dropbox/cloud_work/papers in progress/death_free_market/Supplementary Material/Figures/energy_government_model.R')
source('~/Dropbox/cloud_work/papers in progress/death_free_market/Supplementary Material/Figures/energy_government_model.R')
source('~/Dropbox/cloud_work/papers in progress/death_free_market/Supplementary Material/Figures/energy_government_model.R')
View(sim_best)
source('~/Dropbox/cloud_work/papers in progress/death_free_market/Supplementary Material/Figures/energy_government_model.R')
library(magrittr)
library(here)
library(ggplot2)
library(gridExtra)
library(data.table)
library(ggrepel)
text.size = 10
dir = here()
setwd(dir)
source("functions.R")
# best-fit model
############################################################################
# best fit hierarchy model
wd = gsub("Figures", "Hierarchy Model/results", dir)
setwd(wd)
sim_best = fread("government_model_fit.csv")
names(sim_best) = c("energy_pc", "government_frac")
sim_best$government_frac = sim_best$government_frac*100
# best fit parameter
gov_nfirms_best = fread("government_nfirms_best.csv")
gov_label = paste("Model government\n= top", gov_nfirms_best, "firms")
# model loess
loess_mod = bin_log_loess(sim_best$energy_pc,
sim_best$government_frac,
adjust = 0.1)
log_loess(sim_best$energy_pc, sim_best$government_frac)
# model loess
loess_mod = bin_log_loess(sim_best$energy_pc,
sim_best$government_frac,
adjust = 0.1)
source('~/Dropbox/cloud_work/papers in progress/death_free_market/Supplementary Material/Figures/energy_government_model.R')
source('~/Dropbox/cloud_work/papers in progress/death_free_market/Supplementary Material/Figures/energy_government_model.R')
hist(sim_best$energy_pc)
hist(log(10sim_best$energy_pc))
hist(log10(sim_best$energy_pc))
hist(log10(sim_best$energy_pc), breaks = 100)
summarise_coord(sim_best$energy_pc)
summary(sim_best$energy_pc)
